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Abstract. We study the asymptotic behaviour of the heat content on a com- 
pact Riemannian manifold with boundary and with singular specific heat and 
singular initial temperature distributions. Assuming the existence of a com- 
plete asymptotic series we determine the first three terms in that series. In 
addition to the general setting, the interval is studied in detail. 
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Dedication 

This paper is dedicated to Stewart Dowker. Stuart has made many contributions 
to the field of spectral geometry and remains active in this area [2J [T3] . The second 
author has been honored to have been a collaborator with Stewart [16l [17] . We 
hope this paper serves as a fitting tribute to our colleague and friend. 

1. Introduction 

Let M be a compact Riemannian manifold with smooth boundary dM, and let 
5 denote the geodesic distance to the boundary. Let and ip2 be smooth functions 
on the interior of M. Then ip\ will represent the initial temperature of M and ip2 
will represent the specific heat of M. Since M is compact and 8M is smooth the 
distance function is smooth near dM. We have to assume that S ai ipi and 5 a2 ip2 
are smooth on a closed collared neighbourhood of dM. The parameters a\ and a.2 
control the growth or decay of ipi and ip2 near dM. Let D be an operator of Laplace 
type on M. Impose Dirichlet boundary conditions to define the realization of D. 
Let dx be the Riemannian measure on M. Since dM is smooth, the corresponding 
Dirichlet heat kernel pu{xi, #2; t), x% £ M, X2 G M,t > vanishes linearly in 
S(x%) and 8{x2) near dM. We suppose a\ < 2 and a.2 < 2 to ensure convergence 
subsequently. Let e~ tD be the fundamental solution of the heat equation for the 
Dirichlet Laplacian. Then 

ui(-,t) ;= e~ tD ipi = \ p M (;x;t)ipi(x)dx 

JM 

represents the temperature of the manifold for t > 0. The heat content Q is defined 
by 

(2(^1,-02, £>)(*) : = / ui(x;t) ■ ip2(x)dx 

JM 

I PM(xi,X2;t)ip 1 (x 1 )ip 2 (x2)dx 1 dx2 ■ 

M JM 

Conjecture 1. Let a\ + a.2 ^ a,\ < 2, 02 < 2. There is a complete asymptotic 
series as t ], 

00 00 

n=0 3=0 
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where the /3„ , n = 0, 1, • • • are regularized integrals of local invariants over M and 
where the (3f M , j = 0, 1, • • • are integrals of local invariants over the boundary. 

Remark 1. It is convenient to let ot\ and cti be complex as we may then use 
analytic continuation. For 3?(ai) << and SR(«2) << 0, 

™ ! hi 

The values of /3„ for more general values of a\ and a 2 may then be obtained as 
regularized integrals as discussed in [TO]. We omit the technical details concerning 
the requisite regularizations in the interests of brevity as they will play no role in 
our analysis. 

The heat content has obvious physical relevance and the invariants fif M , which 
reflect the asymptotic behaviour ast 10, relate the geometry of M to the underlying 
physical properties of M. Much of the previous work in the field has been devoted 
to the computation of the invariants /3^ M in the smooth setting (ai = 0, a 2 = 0). 
They were originally studied for the scalar Laplacian with ipi = fa = 1 [21 HI CO] ■ 
Subsequently, general initial temperatures and specific heats were investigated - 
see [5J |TU1 [T31 [TH1 [201 (HI 122] and the references contained therein. Other boundary 
conditions (Neumann, Zaremba, etc.) have been considered [31(5]. The growth of 
the coefficients /3^ M has also been of interest [JJ [71 [23] ~ see also [12] for related 
work on the heat trace asymptotics. The case where ipi is singular («i > 0) but ip2 
is smooth (a 2 = 0) was studied previously [51 UOj. The current paper is devoted to 
the study of the invariants (3j M in the doubly singular case. 

The special case of a ball of radius a in E 3 is well understood. The following 
result was proved in [6]. 

Theorem 1. Let B„ = {i £ I 3 : \x\ < a}, and let D be the Dirichlet Laplacian 
acting in L 2 (B a ). Lf ct\ < 2, a.2 < 2, ct\ + cv.2 > 3, J G N then there exist coefficients 
bo, b\, ■ ■ ■ depending on a\, a.2 only such that for t \. 

Q{5- a \8- a \D){t) 

— Aire r7 2 /( 1 - Q i- Q 2)/2 _ 4_c„ i A- r ,)nf( 2 - a i~ a 2)/ 2 

J \ L ) 

+47r CQl _ 1 , Q2 _ 1 i( 3 - Q i-^)/ 2 + J2 6 J a 3 - J '- ai - aa * I '/2 + 0(t< J+1 >/ 2 ), 

where 

c aua2 = 2- Q i-^7r- 1 /2r((2 - ai - a 2 )/2) 

/- 1 (2) 
x / (p- ai +p- a2 )({l-p) ai+a2 - 2 -(l + p) ai+a2 - 2 )dp, 
Jo 

and 

b = -87r((ai + a 2 - l)(ai + a 2 - 2)(«i + a 2 - 3)) _1 , 61 = 0, 
b 2 = 87raia 2 ((ai + a 2 + l)(«i + a2)(ai + a 2 - I)) -1 , b 3 = 0. 

It is convenient to use a standard formalism to describe the invariants (3j M in the 
general setting. Let D be an operator of Laplace type acting on the space of smooth 
sections to some vector bundle V over a Riemannian manifold (M,g). Choose a 
local system of coordinates (x 1 , ...,x m ) for M and a local frame for V. We adopt 
the Einstein convention and sum over repeated indices. Let ds 2 — g^dx^ o dx v 
define the Riemannian metric and let g^ v be the inverse matrix where 1 < /a, v < m. 
We may then express: 

D = -(g^Idd Xfl d Xu +A^d IM + B) 
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for suitably chosen endomorphisms A^ L and B oi V . If V is a connection on V, 
we use V and the Levi-Civita connection to covariantly differentiate tensors of all 
types and let ';' denote multiple covariant differentiation. If tpi is a section to 
V which is smooth on int(M), let Vi;/*" ^ e the components of V 2 ^i. If E is an 
auxiliary endomorphism of V, we define the associated modified Bochner Laplacian 
by setting: 

D(g,V,E)^ :=-ffer^i' 
Let T^ a and r /J „ <T be the Christoffel symbols of the Levi-Civita connection. Then 
(see, for example, the discussion in [S]): 

Lemma 2. If D is an operator of Laplace type, then there exists a unique connec- 
tion V on V and a unique endomorphism E of V so that D = D(g,\7 , E) . The 
connection I- form uj of V and the endomorphism E are given by 

uj^ = ±{g^A» +g° £ T aetl Id), 

E = B - g^ u (d Xiy u^ + cj m cj„ - uj a T f _ lly a ) . 

The specific heat ip2 is a section to the dual vector bundle V. We use the dual 
connection on V to covariantly differentiate ip2- Note that the connection 1 form 
(i„ for V is the dual of — lo u . Thus 

Va^ - d x „ - \{ 9ilv A» + g^T^ld) . 

Near the boundary, choose an orthonormal frame {ei,...,e m } for the tangent 
bundle of M so that e m is the inward unit geodesic normal. Let indices a, b range 
from 1 to m — 1 and index the induced orthonormal frame {ei, e m _i} for the 
tangent bundle of the boundary. We let ':' denote the components of tangential 
covariant differentiation defined by V and the Levi-Civita connection of the bound- 
ary. Let L a b '■— g(^e a eb, e m ) = ^ abm be the components of the second fundamental 
form. The difference between ';' and ':' is then measured by L. Let D be the dual 
operator of Laplace type on V. The following relations will be useful subsequently: 

Dipx = -(il>i; aa + ipi ;mm - L aa 1pi- m + Eipi), 

Dlp2 = -{^2:aa + ^2;mm ~ L aa 1p 2 - m + Elp 2 ) ■ 

We expand ipi and ^2 near the boundary of M in the form: 

oo oo 

where S7 em ip{ = an d ^e m ^2 = 0- We shall usually be working with scalar 
operators and can choose local sections s and s so that V em s = and V em s = 0. 
We may then express ipi = "Jis, ip\ = ^>\s, tp 2 — &2S, and ip^ — where \&j 
and ^2 are smooth functions defined on the boundary so that we have the modified 
Taylor series 

00 

5 a ^ i (y,5)^^2^i(y)6 j for i = 1, 2 . 
3=0 

For the Laplacian, the bundles and connections under consideration are trivial so 
this formalism is unnecessary. However, for more general operators, the connections 
in question are not flat and this formalism is essential. We shall be using the method 
of "universal examples" in what follows. It is a peculiar feature of this method that 
even if we were only interested in the scalar Laplacian for a smooth bounded domain 
in M. m , it would be necessary to deal with quite general operators as we shall see 
presently while proving Lemma [TT] in Section [3] 
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Let Ric be the Ricci tensor of M, let t be the scalar curvature of M, and let 
dy be the Riemannian measure of dM. Section [3] is devoted to the proof of the 
following result: 

Theorem 3. Let ot\ + 012 £ Z, a% < 2, 012 < 2. Assume that Conjecture^ holds. 
Let c ait0 . 2 be as given in Equation 0). Then 

Pq M — Jqm C ai-a 2 Ipl ^d-U, 

+c ai , a2 -ii>l^\}dy, 

02 M = /gjn{ C ai-2,a 2 ^1^2 — \{ c a x -2,a 2 + C ai - l,a 2 -1 )L aa 4>l ^2 

+Ca ua2 Elpi1p2 + C ai , a2 -2i>l1p2 ~ | (Cc*i-l,a 2 -l + C Ql:Q2 _ 2 )iaaV'lV ; 2 
+ ( — jC Qi _2,q 2 — \ c a 1 ,a 2 — 2 + \ c a- s _,a 2 ){Laa.I J bb + RiCmm)^?^ 
-C ai , a2 1pi:aP20:a + 0rip°1p° + C ai -l,a 2 -l^l V>2 

+ (gC ai -2,a2 + \ c otx.a 2 -2 + \c ai -\,a 2 -\ ~ ^C QliQ2 )L aa Lf,b'0l'02} ( ^2/ ■ 

For the ball B a in R 3 , L aa = 2a _1 , L aa Lbb = 4a -2 and L aa Lbb = 2aT 2 . Hence 

/3f Ba = -47ra(c Q , 1 _i !Q , 2 +c aijaa _i), 

and the coefficient of £( 2 - Q i- Q 2)/ 2 agrees with Equation (TTJ) in Theorem|31 Similarly, 
the next term in the series given by /3 2 Ba is consistent with Theorem [3l 

We observe that the V- function in the expression for c ait0 . 2 which is given in 
Equation ^ has simple poles for «i+ct2 £ {2, 4, 6, • • • }. Furthermore the integrand 
with respect to p equals for ct\ + 012 = 2. It is easily seen that this singularity 
is removable. On the other hand the integral with respect to p is finite only for 
ati < 2,ct2 < 2 and a\ + a 2 > 1. This suggests that the j th term (j = 1,2,3) in 
Equation (TTJ) will take a different form for ai + 012 — j- This is indeed the case for 
an interval in R. Let a > 0, and let xi,X2 be non- negative C°° functions on R + 
defined by 

1 if < x < ei, 2 , 

if 2! > £3,4, 

where < £1 < £3 < a/2 and < £2 < £4 < a/2. Wc shall establish the following 
result in Section [2j 



xiA x ) 



Theorem 4. Let u\ < 2, 012 < 2, u\ + U2 = I. If 1 10, then: 

P[o.a](xi,x2;t)xi(S(xi))x2(S(x2))S(xiy ai S(x2y a2 dx 1 dx2 

[0,a] 2 

= log(£ 2 /t) +7 + 4 log(2 1 / 2 - I) + 4 log 2 

+2/ Xl(x)x2(x)x- 1 dx+ dqq-^l + q 2 ) .* 

J[£,a/2] .7 [0,1] W 

xjca+^Ai-g^ + ai-ffva+gjr 

-2(I-a)(I + a 2 )- 1 / 2 | + 0( i i/2 logi ) 7 

where P[o, a ] %2', t), x\ € [0,a],X2 G [0,a],i > is £/ie Dirichlet heat kernel for the 
interval [0, a], 7 is Euler's constant, and e = min{£i,£2}. 

We note that the £-dependence in the right hand side of Equation (J3J is fictitious. 
Since xi( x ) = X2(x) = 1 for < x < e, we have that 

log(£ 2 /i) + 2/ xi{x)x2(x)x~ 1 dx = 2 / Xi{x)X2{x)x~ x dx, 

J[e,a/2] J[Vt,a/2] 
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which independent of e for < t < e 2 . We also note that the leading term in 
TheoremjUjibes with Theorem 1.4 (2) in [TU] since the volume of the volume of the 
boundary of the interval [0, a] is equal to 2. This supports the following. 

Conjecture 2. Let M be a compact Riemannian manifold with smooth boundary 
dM, and let S denote the distance to the boundary. Let ai < 2,a2 < 2,ai+a 2 = 1; 
and let \i an d Xi be smooth functions on R + with support contained in an interval 
[0,6], and equal to 1 in a neighbourhood o/O, and where b is such that S is smooth 
on the collar dM x [0,6]. // 1 J, then 

Q{5-^xio5,5-^ X 2o5,D)(t) = 2- 1 f dylogt + o(logt). 

JdM 

We note that Theorem @] and Conjecture [H include the cases where either 1 < 
ai < 2 or 1 < a 2 < 2. This requires more care in the proof of Theorem 2] than the 
case where both a\ < 1 and a 2 < 1- 

2. The proof of Theorem [J] 

The first step in the proof of Theorem [4] is to reduce the calculation on the 
interval [0, a) to a calculation on the half-line K + = [0, oo). We have the following: 

Lemma 5. Let a\ < 2, a 2 < 2, a.\ + cli = 1. // 1 J, then 

P[o, a ] (xi 1 x 2 ;t)xi{S{x 1 ))x2(5(x 2 ))S{x 1 y ai S(x2y a2 dxidx 2 

PK.+ (xi,X2;t)xi(x 1 )x2{x2)xi ai x 2 ~ a2 dxidx 2 ( 4 ) 

2 

+ 

+ ( e -(l-'))- 2 /(4t)) j 

where n = a — £3 — £4 and 77 = max{ai/2, 0:2/2}. 

Proo/.Without loss of generality we may assume that ai > ct2- We partition the 
region of integration [0,a] 2 = uf =1 A;, where 

M = [0, e 3 ] x [0, e 4 ], ^2 = [0, e 3 ] x [a - e 4 , a], 

A3 = [a - e 3 , a] x [0, e 4 ], A4 = [a - e 3 , a] x [a - e 4 , a], 
A 5 = A\(uf =1 A i ). 

The integrand in the left hand side of Equation Q is identically equal to on A5, 
and this set does not contribute to the integral. Since 

P[o,a] Oil x 2 ; t) = pp,a] (a - xi , a - x 2 ; t), and 

P[0, a ](xi,a- x 2 ;t) =P[ ,a](a - x u x 2 ;t), 

the contributions of A\ and A 2 to the integral in the left hand side of Equation @ 
are equal to the contributions of A 4 and A3 respectively. Since \x\ — x 2 \ > k for 
(xi,x 2 ) € A 2 , we have by monotonicity of the Dirichlet heat kernel that 



P[o,a]{xi,x 2 ;t) < PR+(x!,x 2 ;t) 

= (AlTt)- 1 / 2 ^ e -( K i- a; 2) 2 /(4t) _ e -(x 1 +x 2 ) 2 /(4t) 
= (4 7 ri)- 1 /2 e -(a:i-a;2) 2 /(4t) (l _ e -Xi<ts a /t\ 

< t- 3 / 2 x lX2 e- K2 /( 4t l 



(5) 



Hence the contribution from A2 to the integral in the left hand side of Equation (j4|) 
is bounded from above by 
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The contribution from A\ to the integral in the left hand side of Equation ^ is 
bounded from above by 

5+ (xi , x 2 ; t)x\ (xi )x2 (x 2 )x^ ai x^ 012 dxidx 2 

p R + (xi , x 2 ; t)xi (xi )x2 {x 2 )xi ai x^" 2 dxi dx 2 . 



This completes the proof of the upper bound. 
To establish the lower bound we note that 



P\o,a}(xi,x 2 ]t)xi(S(xi))x2(S(x 2 ))S(xi) ai 6(x 2 ) a2 dxidx 2 

[0,q] 2 



~ 2 J J p ^' a ^ Xl,X2 '^ Xl ^ Xl ^ X2 ^ X2)Xl aix 2 a2 dxidx 2 . 

It is well known that the Dirichlet heat kernel for an open set fi G R m has the 
following probabilistic representation. 

pn{x\,x 2 ;t) = p Rm (xi 1 x 2 ;t)Piob Xl . X2 [B(s) G f2,0 < s < t], 

where (B(s),0 < s < t) is a Brownian bridge on R m . For ft — [0,a] G R and for 
x G [0,a],j/S [0, a] we have that 

P[o,a]{xi,x 2 ;t) = p R (xi,x 2 ;t)Prob 2;i , a:2 [0 < B(s) < a, < s < t] 
= Pr{%i, x 2 \ t)(Prob ;EliX2 [0 < B(s),0 < s < t] (6) 
-Prob Xl ^ 2 [(0 < B(s), < s < t) A (max < s < 4 B(s) > a)}). 
By Holder's inequality we have for r\ G (0, 1) 

PR(xi,x 2 ;£)Prob 

xi,x 2 [(0 < B(s), < s < t) A (maxo< s <t B(s) > a)] 
<»(xi,x 2 ;i)(Prob Xl!:E2 [0 < S(s),0 < s < i])'' ^ 

x(l-Prob xl , !C2 [S(s) <a,0<s<t]) 1 -" 
= (p R +(xi,x 2 ;^)) , '(pR(xi,x2;t) -p(-oo,a](zi,X2;i)) 1 ~' ? . 

Since 

P(-ao,a](«i^a;*) = (4^r V2 (e- (ai - s 2)2/(4t) - e -(**-*i-*»>VC«)) , (8) 
we have by Equation ([5]), Equation (JTJ), and Equation © that 

P[0,a](zi,X 2 ;i) > PR+(xi,X2',t) 

-{p^{ Xl ,x 2 -t))^AT:t)-^l 2 e-^-^ 2a -^- x 
By the last inequality in Equation ((5J) 

(p R+ (x 1; x 2 ; i))^^)-^-")/^-^-")^-^-^) 2 /^) 

<t-^{x 1 x 2 ye- ( - 1 '^ 2a -^- X2 ^' l - it l 

Integrating the above right hand side with respect to Xi(xi)x 2 (x 2 )x||~ ai x 2 ~ Q2 dxi(ix 2 
yields a bound 

/ xi(^l)^i" ai ^i / Xi(^)a;r Q2 ^2e^ (1 - r ' )(2a - a;i - a;2)2/(4t) 

•/[0,e 3 ] J[0,e 4 ] 

< t -i-T, e -(l-^)(2 -e 3 -e 4 ) 2 /(4f) /" Xl ( Xl ) x ^ d Xl f X l{x 2 )x\~ a2 dx 2 

J[0,e 3 ] J [0,64] 

= 0( e -( 1 - r ')( Q - e 3- e '4) 2 /( 4t )). 
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Note that since 2 > a± and 77 = Qfi/2, x^ " 1 = 2^ ai ^ 2 is integrable at 0. Since 
1 = a>i + a 2 < 2ai we have that u\ > 1/2 > 0. Hence x!p Q2 = x 2 1+ ( 3ai / 2 ) j s a j so 
integrable at 0. This completes the proof of the lower bound. □ 
In order to prove Theorem 0] it clearly suffices to prove the following. 

Lemma 6. Let a 2 < ct\ < 2, a,\ + a 2 = 1. If 1 1 then 

p R + (x\ , x 2 ; t)xi {xi)x2 {x 2 )x^ ai x^" 2 dxidx 2 

2- 1 log(e 2 /i) + 2" 1 7 + 2 log(2 1/2 - 1) + 2 log 2 

+ / Xi{x)X2{x)x- 1 dx + 2- 1 ( dqq-^l + q 2 ) (9) 

J[e,a/2] .7 [0,1] 



-2(1 - g)(l + q 2 y 1/2 ^j + 0(t x ' 2 logt). 



Proof. Define 



C = {{xi,x 2 ) eM. 2 + :x 2 1 +xl>e 2 ,0<x 1 < e 3 ,0 < £2 < e 4 }, 
Ci = {(xi,x 2 ) € C : |xi - x 2 \ < a}, 

where a G (0, e/5) will be chosen later on. The left hand side of Equation (O can 
be written as B\ + B 2 , where 

By = I I p R +(x 1 ,x 2 ;t)x^ ai X2 a2 dx 1 dx 2 , (10) 

J ,/K2_n{2:f+a;|<£ 2 } 

d k + (xx , x 2 ; t)x\ (xi)x 2 (x 2 )xi ai X2~ a2 dxidx 2 . 

c 

To estimate B 2 we first consider the contribution from the set C\C\. We have 
by Equation ([5]) that 

P &+ (x 1 ,x 2 ;t) < t-^ 2 x lX2 e-^- x ^l^ < t- 3 ' 2 x lX2 e-^'^\ (x u x 2 ) € C\C V 
Consequently, 

-{x x ,x 2 ;t)xi{x 1 )x2ix 2 )x^x- 2 a2 dx x dx 2 < Kt- 3 / 2 e- a2 /^ (11) 



'C\Ci 
where 

K = J J Xi(xi)X2(x2)x{~ ai xl~ a2 dxidx 2 . 

On C C\{\xi — x 2 \ < e/5} we have that x 2 — > X2{x 2 )x 2 ~ a2 is C°°. Hence there exists 
L depending on e,a 2 and on X2 such that \x2(x2)x2~ a2 — X2(i)£;T Q2 | — L\x\ — x 2 \. 
It is easily seen that both x\ > e/2 and £2 > e/2 on C n {|xi — seal < e /5}- Since 
the Dirichlet heat kernel on R+ is bounded from above by i -1 / 2 we have that 

p R+ (xi,a;2;i)xi(a;i)a;j; Ql L|a;i - x 2 \dxidx 2 

r r (12) 
< L(2/e) a H- 1 / 2 / |xi -x 2 |rfa;irfa;2 < 2aL(2/e) Ql r 1/2 cr 2 . 



We now choose a 2 as to minimize t 3 / 2 e ct2 /( 44 ) + t 1 ^ 2 a 2 , i.e. 

cr 2 = 4ilog(i- 2 ). 
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This gives that for t sufficiently small the right hand sides of Equation (TTTT) and 
Equation (TT2")) are O^ 1 / 2 ) and 0(t Y l 2 log(i -1 )) respectively. We conclude that 

B 2 = J J p R +(xi,x 2 ;t)xi(xi)x2(xi)x^ 1 dxidx2 + 0(i 1/2 log^ 1 )). (13) 

We now write 

C x = (Ci n {x\ > e 2 /2}) U dinixl < e 2 /2}) = C 2 UC 3 . 

Since x\ > e/2 on Ci we have that the integrand in the first term in the right hand 
side of Equation (IT31) is bounded by 2e~ 1 t~ 1 / 2 . Hence 

Ps+(xi, x 2 ;t)xi(xi)x2{xi)x 1 ~ 1 dxidx 2 < 2e~H~ 1 < 2 \C 3 \, 

c 3 

where | ■ | denotes Lebesgue measure. It is easily seen that | C3 1 < er 2 /2. Conse- 
quently, 

0<JJ Pm,+ {xi 1 x 2 ;t)xi{xi)x2{xi)x^ 1 dxidx 2 <e~ 1 t~ 1/2 o 2 , 

and so the contribution from C3 to the integral in Equation ([T3"| is 0{t x l 2 log(t -1 )). 
Furthermore by monotonicity of the Dirichlet heat kernel 



Hence 



p K +(xi,x 2 ;i) < pm{xi,x 2 ;t) . 
Pr+ {xi,x 2 ;t)xi(xi)x2(xi)x^ 1 dx 1 dx 2 



c 2 

< / / PR(xi,x 2 ;t)xi(x 1 )x2(xi)x 1 ~ 1 dx 1 dx 2 

J J{x 2 >e 2 /2} 

Xi(xi)x2(xi)x^ 1 dx 1 . 

/V2,a/2] 

To obtain a lower bound for the contribution from C2 to the integral in Equa- 
tion CGD we first observe that (4 7r t)-i/2 e -(^i+^2) 2 /(«) < t -i/2 e ~e 2 /(4t) and x > ( j 2 

for (£1,0:2) G C 2 . Therefore 

< / / (inty^e-^+^/^xiix^ix^dx^ 
J Jc 2 

< 2e-H- 1 / 2 e^ 2 /^\C 2 \ < 2a 2 e- 1 t- 1 / 2 e~^/ { ^ 



Finally 



c 2 



{xi,x 2 ;t)xi{xi)x2{xi)x~ 1 dxidx2 



> I p R (x 1 ,x 2 ;t)xi(x 1 )x2(xi)x 1 1 dxidx 2 

'{x 2 >e 2 /2} 

Pr{xi, x 2 ;t)xi(xi)x2(xi)x 1 ~ 1 dxxdx 2 

{\x 1 -x 2 \>(j}n{x 2 1 >e 2 /2} 

Xi(xi)X2(xi)x^ 1 dx 1 

/V2,a/2] 

Pr(xi,x 2 ; t)xi(xi)x2(xi)x 1 ~ 1 dx 1 dx 2 . 

{\x 1 -x 2 \>a}n{x 2 >e 2 /2} 
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Moreover 

Pm (x 1 ,x 2 ;t)dx 2 < f {Ai:t)- 1 ' 2 e-^- x ^'^dx2 

{\x 1 ~x 2 \><y} J {\xi-x 2 \>(t} 

= 47T- 1 / 2 i 1 /2 (J -l e ^ 2 /(4t) = 0(t 2 ). 

Putting all this together gives that 

B-2 = I xi(x) X 2(x)x- 1 dx + 0(t 1 / 2 log(t- 1 )) 

J[e/\/2,a/2] 

Xi(x) X 2(x)x- 1 dx + 2- 1 log 2 + 0(i 1/2 log^ 1 )), 

,o/2] 

since xi{ x )X2{ x )x^ 1 = x^ 1 for < x < e. 

In order to obtain the asymptotic behaviour of B\ in Equation (110|) , we introduce 
polar coordinates x = {At) 1 / 2 p cos 6, y = {At) 1 / 2 p sin 9 to find that 



B x = n- 1/2 [ d9{cos8)- a {sin9) a -' 



X / ^ (e -^(l-si n (20))_ e -^(l+sin(20) ):) _ 
J[0,e/(4t)V2] 

A further change of variable 9 = <j) + ir/A yields that 

— 1 /'/-.^o A civ, ±\a — l ' 



= 1/2 f /(cos^ + sin^ (coB^-sin^\ 

7[o l7 r/4] V (cos0-sin0)« (cos + sin 0)" J 

X [ dp ( e -2p 2 (-n0) 2 _ e -2p3(cos«n = B3 + B4 + Bs 

./[O.e/Mt') 1 / 2 ! V ' 



where 



B 3 = (2/tt) 1 / 2 / # ( ^ Z D111 ^ + (C ° S ~ Sln 0) ° 1 - 2 
/[0,tt/4] V (cos 0- sin 0)" (cos</> + sm(/)) Q 

. / dp(e- 2 P 2 ^ n ^ 2 _ e -2p 2 (cos0) : 

'[0,00 ) ^ 

2" 1 / d0(cos0) _1 (sin0) _1 

J[0,ir/4] 

'(cost/) + sine/))" -1 (cos — sin q. 



{cosip — sm<p) a ~ 1 (cos0 + sin0) Q 

dqq-^l + q 2 ) X ( ((1 + q)/{l - q))^ 1 



2(cos<^-sm</>) (14) 



[04] 

+ ((1 -<?)/(! + <z)r- 2(1 -gXl + g 2 )- 1 / 2 



B 4 = -(2/tt) 1 / 2 



[0,tt/4] 



(cos^ + sin(/)) Q 1 (cos — sin 1 
^ (cos <^> — sin </>)" (cos</> + sin0) Q 



x 



dp f -2p 2 ( S i„0) 2 _ e -2p 2 (cos0) 2 \ _ 
/(4t) 1 /2, oc ) V / 



and 





4 / ^( 


/[0,7r/4] 


i[0,e/(4t)i/2] \ 



-2p 2 (sin0) 2 — 2p 2 (cos 



e -2 P (COS?, ^ 
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We have used the standard change of variables tan <f) = q to obtain the last identity 
in Equation (TT4"| . 

In order to find the asymptotic behaviour of B§ as t \, we first consider the 
contribution of the second term in the integrand with respect to p in Equation (|15p . 
and write 

-(8/tt) 1 / 2 / # / dpe- 2 P 2( ~ cos ^ 

J[0,tt/4] J[0,e/(4t)i/2] 

d^cos^)- 1 + {8/it) 1 / 2 f d(j> ( dpe- 2p2(cos ^ 2 

[0,7r/4] J[0,K/i] J[<L/{ity/ 2 ,Oo) 

= log(2 1 /2_i) + ( e - 2 /(5t) ) . 

The contribution of the first term in the integrand with respect to p in Equation f|15p 
is calculated as follows: 



(sm t; 

,,2 .2 



( 8 /-) 1/2 /[0,V4]^/[0,e/(4 t )V 2] ^e- 2 P ( S 
= (8/^) 1/2 / [(W 4] # /[0,e/(4 t )V 2] rfpe- 2 " 2 ^ ^ 

+W]^(( sin ^)~ 1 -^ 1 ) ' ' 

The third term in the right hand side of Equation (TlT)]) is 0(e~ £ /( 5t )). The second 
term in the right hand side of Equation $W\i is equal to log(2 1 / 2 — 1) + 3 log 2 — log n. 
The first term in the right hand side of Equation (|16[) equals 



i[0,7T£/(32t)l/2] J[O,0] 

(4/ir) 1/2 (log(7re/(32i) 1/2 ) / dpe"" 2 

V 7 i[0,7T£/(32t)l/2] 



-(4/tt) 1 / 2 / ^(log^e- " 
= 2- 1 \og(e 2 /t) + logTr + 2 -1 7 - 3 • 2" 1 log 2 + 0( e - £2 /^), 

where we have used Equation (4.333) in [IS] together with 

/ #(log^)e- 02 = / (^(log^e'^ +0(e" £2/(5t) ). 

J[0,7T£/(32t) 1 /2] J[0,00) 

We find that 

B 5 = 2- 1 log(e 2 /i) + 2" 1 7 + 21og(2 1 / 2 - 1) + 3 • 2" 1 log 2 + 0(e~ £2/(5t) ). 

In order to estimate B4 we first note that by expanding sin <fi and cos <f> around 
we have that 

(cos^ + sin^)"" 1 (cos0 - sin (ft)"" 1 _ ^ = Q 2 
(cos </> — sin (cos</> + suk/>) q 
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Furthermore for cf> € [0,7r/4], 



Hence 



e 2p 2 (sin0) 2 _ e -2p 2 (cos< 



e/(4t) 1 /2, 00 ) 

< (4i) 1/2 e- 1 



/(4t) 1 / 2 ,oo) 



dp 



-2p 2 (sin (; 



J[e/(4f) 1 /2 s0O ) V 



-2p 2 (cos 



-2p 2 (cos <: 



1/2-1 



[0,00) 

-2 



e -2p-(sin0) 2 _ e -2p 2 (cos<; 



i 1 /V 1 ((sin0)- 2 -( 



COS < 



\B 4 \ < Tf-^fW 1 / #(sin0)- 2 (c 



COS0 + S11K 



+ 



cos0 — sine 



— - 2{{cos <f>y - (sin^) 2 ) 



(cos<p — sin^)"" 1 (cos + sin < 

We see that the integral with respect to <f> converges both at <fi = and at (j> = ir/4. 
We conclude that B 4 = C^i 1 / 2 ). □ 



3. The proof of Theorem [3] 

We shall assume that 5ft(ar) << and 3?(a2) << and then apply analytic 
continuation to establish the general case. We shall also assume that a% + £ Z 
to ensure that the interior and the boundary terms do not interact. The invariants 
are given by local formula. Standard arguments using dimensional analysis 
yield the following result; as these arguments are by now standard (see, for example, 
the discussion in [5]), we omit details in the interests of brevity. 

Lemma 7. There exist universal constants e\ „ so that: 

P 9 2 M = / 9 M{4 1 ,a 2 (V'?,V' 2 )+4 1 ,a 2 (Wl 1 ^ 2 )+4 1 ,a 2 <^l,V' 2 ) +< l0a <^> 



+4!,Ji«Wl° ) €> +4 2 1 ,a 2 (^:a,€ :a ) + 4U<^1,V> 2 °> 

+4 4 1 , Q2 ^i 1 > 2 1 )}dy- 



Remark 2. We note that e„ 



is given by Equation ©. 



There is a basic symmetry which is useful. Let e~ tD denote the fundamental 
solution of the Dirichlet Laplacian and let D be the dual operator on the dual 
vector bundle V . The lemma below follows immediately from the identity 



Q(i/>i,i>2,D)(t) 



{i> 1 ,e- tu i )2 )dx = Q{il>2,il>i,D){t). 



IM JM 

Lemma 8. Adopt the notation of Lemma^ 



F o _ F o F i 

cti,ct2 a 2 ,ai' Qi,tt2 

~5 _ p8 F 6 

-.11 _ _11 p 12 

ai,a 2 a 2 ,ai' t ai,a2 



£ 3 £ 2 

a 2 ,ai ' ai,a 2 

e 6 e 9 

rl2 



"Ot2 -Ql ' 



2 4 7 

F f = F 

Q2.Ql' a'l ,ft2 Q2.Ctl' 

r 9 ' p \6 -10 ' 

t Q2,Ql! CUl,CK2 Q2,Q!l' 

13 _ -13 ^14 _ ^14 



ai,Q2 



"02 ,Q!-1 ' 



Next, we consider some product formulae: 
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Lemma 9. Suppose that M = M\ x M 2 , that gu — 9m x + 9m 2 , that dM\ = 0, and 
that Dm = Du x + Dm 2 where Dm ± and Dm 2 are scalar operators of Laplace type 
on Mi and on M 2 , respectively. Suppose that ip M = V'f^f 2 and ip^ = ^2' ^2 2 
decompose similarly. Then 

(a) ^,^,D M )(t)=^\^\D M J(t)-p(^,^,D M2 )(t). 

(b) f BM « >aa «, ^,D M )dy = £ 2 „ +j=fc J Ml 1 , (A* )"^ 2 Ml )dx Ml 

(c) TTie universal constants e l ai CC2 are dimension free. 

(d) e e - e° e 13 - and e 12 - -e° 

Proof. Assertion (a) follows from the identity e~ tL>M = e~ tDM ~t eT tD ^ 2 and Asser- 
tion (b) follows from Assertion (a). If we take M\ = S 1 , Dm x — —dg, tp^ 11 — 1, and 
ip% h = 1, we have that fi{ij) Ml ,D Ml ){t) = 2ir. This then yields the identity 

5 3M /■»/,M 2 „/,M 2 n \^, - / odM 2 /„/,M 2 j,M: 



dM JdM 2 

Assertion (c) now follows. We take M 2 = [0, 1] and D 2 = —d^. We take 

ijjMi = ^2 = near r = 1, 

^m 2 = r - a2 and ^m 2 = r - ai near r = Q 

Since the structures on M 2 are flat, we have = = for A: > while 

v> 2 = ^ 



I at r = l 1 
\ 1 at r = J ' 



Consequently, 

{0 if r = 1 and k > 

if r = and k>0 

e° aua2 if r = and k = 

As the second fundamental form vanishes, the distinction between ';' and ':' dis- 
appears, and we have Diip^f 1 = —(tp^aa + Efa)- Calculating on the interior then 
implies that 



mr^2\D Ml ) = / (^,^ a + E^)d Xl . 

We may therefore use Assertion (b) to derive the following identity from which 
Assertion (e) will follow: 

/ $™ ua2 {^,i>¥,D M )dy = el ua2 f (^\^L + E^)d Xl . □ 

JdM J Mi 

We continue our study by index shifting: 
Lemma 10. 

F \ _ F F 4 _ _5 _ _2 

c ai,ct 2 c qi-1,q 2 i °ai,a 2 °ai— 2.a 2 7 °ai.a 2 c ai-l,a 2 i 

P 3 — P° F 7 — F° F S = P 2 

°a-i,a 2 c ai,a 2 — I? c ai,a 2 c ai,a 2 — 2' °ai.a 2 c qi,q 2 - 1' 

-14 = F 

c ai,a 2 ai — l,a 2 — 1" 

Proof. We assume ■01 and ip2 have compact support near the boundary of M. We 

set 0i := (<5™ 1 V')(5 _ai_ ™ 1 and 2 := (<5 ni ?/; 2 )<S~ a2 ~" 2 for n t e N. We compute: 

y^ t (i+k- ni - ai - n2 - a2 )/2 f p ktni+aitna+aa (^ 1 ,^ l ,D)dy 

o JdM 
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We set k — £ + ni + n2 and equate powers of t to see 

Note that i/^ +ni = and Vv™ 2 = "02- The desired result now follows by taking 
(n 1; n 2 ) = (1,0), (0,1), (2,0), (1,1), and (0,2). □ 

Lemma 11. Lei T" 1 ^ 1 denote the torus with periodic parameters (j/i, y m -i) and 
let M := T 1 "- 1 x [0, 1]. Let f a G C°°([0, 1]) /wwe compact support near r = witft 
/o(0) = 0. Lef 0(r) 6 C°°([0, 1]) /icwe compact support near r = imfft = 1 «ear 
r = 0. Let 5 a G M. >Sef 

ds lf = E a e 2ta{r) dy a o dy a + dr o dr, V2 := 0(r)e" /a(r) r -a 2 ; 

:= - Ea z- 2fJr) (dl + S a d ya ) - dl ^ := G(r)r-^ . 

(a) Ifk > 0, then ! BM ft% ua ^wfo,D M )dy = 0. 

(b) -ie 1 - e 2 - ±e 3 =0 

v u / 2 01,02 ai,02 2 ai,aj 
f c ) _If e 6 ,12 N = q 

Cd) -! e 4 +i £ 6 - i e 7 - e 9 =0 

yyj.j 4 ai,a2 2 ai,Q2 4 ai,a2 0:1,0:2 

f e -) I £ 4 ,1 £ 5 , 1 £ 6 , 1 £ 7 , 1 £ 8 ,10 ,1 £ 14 =0 

VW 8 ai,aj 1 2 «i,«2 4 oi,«2 8 01,012 1 2 oi,« 2 1 01,02 4 01,02 

ff) -£ 9 +£ U =0 

W °ai,a2 T °ai,02 

Proof. We use — d 2 on [0,1] and -Dm on M. Since vanishes near r = 1, 
this boundary component plays no role. Let u(r;t) be the solution of the heat 
equation on [0,1] with Dirichlet boundary conditions and initial temperature 
The parameter r is the geodesic distance to the boundary near r = 0. Since the 
problem decouples, u(r;t) is also the solution of the heat equation on M with 
Dirichlet boundary conditions. The Ricmannian measure 

dx = ^/det gijdydr = fa dydr . 

As V>2 = ee-SJ« r -«2 i ^ 2(fa , = Qr- a *dydr. Since vo^T" 1 - 1 ) = (2 7 r) m - 1 , 

Q{i>ui> 2 ,D)(t)= J u(r;t)^ 2 dx = {2n) m - 1 J u(r;t)e(r) r - a2 dr 

= (2n) m - 1 f3(er- ai ,Or- a2 ,-d?)(t). 

The structures are flat on [0,1]. Since vanishes identically near r = 1 and is 
identically 1 near r = 0, only the term /?o is relevant in computing the boundary 
terms; the f3k, ai ,a 2 vanish for k > 1. 

To apply Assertion (a), wc must determine the relevant tensors. We have: 

Tabm = -f'Jabe 2f " , T ah m = -f a e 2 ^6 ab , 

r amb — fa^ab^ ^ a t L am fa^a-bi 

L a b = F a b m \dM = ~fa^ab, 

u a = l-e 2fa 5 a , Co a = -u a = -\e 2fa 5 a: 

U m = ~\ Ea /a J W ™ = = I Ea fa ■ 

Consequently: 

-R-ambm y((^a^m ^m^a)^o; ^m) — T ac T m 5 C^nTab 
= {-(/a) 2 + /a' + 2(/ a ) 2 }e 2 ^a6, 

Ri Cmm --Ea{/a+(/a) 2 }, 

— I V f"-iV A 2 — I V f ' f ' J. i V f f 
2 Z_/a a 4 zl^a a 4 L-ia,b J aJb ' 2 £—ia,b J a J b 
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— IV f " - I V (5 2 + IV f f 

4 Z_/a a 4 Z_/a,b ^ aJ b' 

We compute: 

tf? = 1, 

^ - {V 9r (r^l)}U = {(Or - | Ea /a)(l)}U = ~§ Ea /a, 

= H( v ^) 2 (^^i)}|oa/ = H( 5 4E Q /:) 2 (i)}la A / 

— 1 ft ft — 1 Y^ f" 
~ 8 2—ia,b J aJ b 4i-ioJo' 

V> 2 ° = i, 

$ = {V dr (^)}\dM = {(dr + \ Ea fa)(e- ^ f «)}\8M = ~\ Ea fa, 

^ = |{(V 9r )V 2 }|aM = U(9r + hEafa) 2 (e-^^)}\0M 

— I Y^ f t f t 1 Y^ f " 

8 2—*a,b J a J b 4 Z^a J a ' 

Considering the term E a /a m /^laf yields Assertion (b), considering the term 
E Q <^a m $2% yields Assertion (c), considering the term Ea /" m /-'faf yields 
Assertion (d), considering the term E a b f'afb i n /^fcf yields Assertion (e), and 
considering the term E a (/a) 2 i n 02 cf yields Assertion (f). 

3.1. The Proof of Theorem \3[ We must now simply trace through the logic 
train. We have computed that: 

£ qi,q2 = c «i— i,a 2 and £ QliQ2 = C aijQ . 2 _i, 

£ a 2 ,c(2 = ~2^a!,a 2 ~^^ai,a 2 ) ~ — 2( c ai — l,a2 c a 1 ,a 2 — l)i 

r 4 _ „ n J ,.7 _ r 

t Qi,Q2 — L ai— 2,a 2 aliu t aiiQ2 i-ai,a2- 2) 

£ Ql,Q 2 = C °l> a 2 an( ^ £ Ql,Q2 = C ai-l,a2-l' 

F 12 _ _6 

CKl ,a2 

=•5 



E 11 

4^-ai— 2,a2 ~A^cti,ot2 — 2 ~t~ ^Cai,Q2' 
£ 10 = -|l£ 4 + ie 5 + l£ 6 + l£ 7 + l£ 8 + l£ 14 1 

{g"CQ 1 _2 i a 2 4(^01—2,02 ^Ql — 1,CK2 — 1 ) 4^Qi,Q2 

"^"gCai ..ct2~ 2 4 (^Qi — 1,Q2 — 1 ~^ ^0:1,02—2) "h 4 ^-ai — 1.0:2 — 1 } 
{( — g H" ~^) c ct\ — 2,a 2 (4 4 4") c ai — 1,q 2 — 1 ( — 8 4") C ai,a2 — 2 

4^1,02}- ^ 
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£ 2 


— I.Q2 


= -!( 




-2,a 2 + 


Cq^ — 1 ,CK2 — 


0- 


£ 2 


,a2 — 1 


= -*( 




— l,a 2 — l 


"T ; (j2 — 


2), 


£ 9 


,C*2 


-l£ 4 


,CK2 


4 c ai, 


+ l£ 6 

a 2 1 2 ai 


)CK2 
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